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NOTES ON FISSION DYNAMICS 


M. Mirea* 


Abstract 


The dynamics of the nuclear fission is a complex phenomenon, being not 
yet described adequately from the theoretical point of view. At present, they 
are not models giving a complete description of the richness of the features 
which characterizes this phenomenon. It is the mean reason for which I called 
this paper Notes on Fission Dynamics, being certain that I will not be able to 
make a global description, but only a picture underlining some particularities. 
So, this mini-overview should be considered only a part of the collection of 
articles treating the nuclear physics, published as a special number in the review 
of the Academy of Romanian Scientists, without an exhaustive character. A 
theory treating the nuclear fission is by excellence based on quantum mechanics. 
That is, a theory concerning the interactions between the smallest pieces that 
constitute a many-body nucleus. But, at present it is not possible to perform 
ab-initio calculations to describe the many-body structure of heavy nuclei which 
undergo fission by starting from fundamental interactions. To make the problem 
tractable, the nucleus as a whole are constrained by some collective parameters, 
associated to some collective degree of freedom. The collective variables are 
forced to vary, leading to a scission of the nuclear system. The response of 
the nuclear system to the external forces is given by the nuclear inertia. The 
mean field potential between the nucleons is obtained after a proper average, 
and then used to solve the Schrodinger equation. The treatments presented in 
this article are based on these simplifying concepts. I will give some examples 
of calculations that include the dissipation and the configuration mixing due 
to radial and angular couplings. The importance of the subject is also briefly 
reviewed. 
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1 Introduction 


In 1939, the fission was discovered. It is understand as the a disintegration process 
of a atomic nucleus leading to two fragments of comparable masses. From the exper- 
imental point of view, the process was evidenced by separating the barium from the 
products of the uranium bombarded with a neutron beam by Hahn and Strassman. 
In the same year, Meitner and Frisch gave an explanation for the appearance of 
these products based on the the liquid drop concepts [1, 2]. The fundamental idea 
can be figured as follows: by overcoming a potential barrier, the initial compound 
nucleus is split into two smaller droplets, setting free a large amount of energy. 
Nobody predicted at the time that a neutron can divide a large nucleus into two 
parts, although for nuclei with masses A >120, the fission process into two nearly 
equal fragments is exothermic. From this discovery, the fission became an field of 
research extremely important, in which unsuspected properties of the nucleus were 
revealed. A lot of phenomenological models were invented to explain the strange 
experimental findings, but the underlining features are the expressions of quantum 
mechanics effects. Therefore, the nuclear fission contributed to the development and 
the validation of the quantum mechanics. 


2 Brief history 


In their seminal paper published in 1939 [3], Bohr and Wheeler offered a very com- 
plex theoretical investigation of the fission process. Their picture is a source for all 
theoretical descriptions of this phenomenon obtained up to now. In order to explain 
the fission cross section properties, they assessed that the neutron induced fission 
cross section is statistically proportional to the number of transition states in the 
saddle point of a multidimensional barrier, at a given excitation energy. It was also 
remarked that the fission process is only one possible channel for the disintegration 
of the compound nucleus, being in competition with the neutron emission and the 
the y de-excitation. Also, from nuclear densities of states evaluations they showed 
that the fission cross section is relatively constant up to several MeV. Some neutrons 
are evaporated at larger excitation energies with some probabilities, and the residual 
nucleus undergoes fission, determining a multi-chance fission phenomena. Within 
these essential features, this picture remained unchanged to the present time. 

The fission showed a multitude of intriguing aspects and an impressive develop- 
ment of the physics of the fission followed up to 1953 when Hill and Wheeler provided 
additional key ingredients used to understand theoretically this process [4]. They 
emphasized that the nucleus is an extremely saturated many body system. There- 
fore the potential felt by a nucleon in the interior of the nucleus is approximately 
independent of the positions of the surrounding nucleons. It implies that it is possi- 
ble to consider the nucleus as a whole, characterizing it with a nuclear shape which 
can be deformed. Accordingly, the single particle potential is essentially collective, 
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being managed mainly by its surface boundaries, as a state of the system as a whole. 
These theoretical assumptions lead to a simplification of the problem, which became 
tractable with mean field models. Instead to treat the interactions between all nucle- 
ons, it becomes convenient to characterize the collective state of a nucleus with the 
help of some constraints or some generalized coordinates associated to macroscopic 
degrees of freedom. In this picture, the nucleons move independently in a mean 
field of a microscopic potential managed by these degrees of freedom. That allows 
us to construct mean fields by using phenomenological prescriptions [5, 6, 7, 8] or 
within self-consistent models [9], able to describe properties observed from experi- 
ments. The constraints associated to the degrees of freedom are allowed to vary in 
time leading to the rupture of the initial compound nucleus and leaving to separated 
fragments. This is essentially our understanding about the fission process. 


In neutron induced fission, a direct reaction between the incident bombarding 
particle and the target can lead to elastic or inelastic scattering, as well as a com- 
pound nuclear system. The compound nuclear system disintegrates after a mean 
life. The simplest treatment is to consider that in the low energy domain, dur- 
ing this short period of time the kinetic energy of the incoming particle is redis- 
tributed among the entire nuclear system. For intermediate or high energies, a 
pre-equilibrium stage should be envisaged, to take into account inelastic processes. 
The compound nucleus formed when a nucleon hits the internal region of the tar- 
get nucleus and the cross section of the process is estimated as a solution of the 
Schrodinger equation for optical model potentials [10, 11, 12]. Once the nucleus is 
considered as thermalized, the disintegration can be treated within statistical argu- 
ments, by calculating the energy widths of the all decay channels. The energy in the 
fission degree of freedom is always smaller or at most equal to the total excitation 
of the compound nucleus. For fixed value of total excitation energy, the fission cross 
section is obtained by summing over the totality of energies in the fission degree 
of freedom. In other words, a sum of all fission channels is performed within an 
integration in which the weight is represented by the nuclear density. 


In the investigation of the fission process, the more impressive discoveries were 
empirical in the early stages. Strange experimental findings at the time contributed 
to the developments of the quantum mechanics, of new concepts in the field. To our 
days, an exhaustive theoretical approach including the richness and the complexity of 
the data behavior issued from the experiments is not available. The nuclear structure 
is still not understood adequately. For example, it is still difficult to reproduce the 
nuclear masses on the entire nuclear table. 


As an example of a strange experimental behavior is the discovery of an isomer 
in americium [13, 14] together with the presence of intermediate structures in the 
fission cross sections [15, 16, 17]. These strange behavior were not explained until 
the postulation of a double humped fission barrier. A such shape of the barrier is 
impossible to be obtained within the liquid drop model, where a constant distribu- 
tion of the nuclear density inside the nucleus should be assumed. The deformation 
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energy should depend on other quantities, as the internal structure of a nucleus. 
The distortion of the nuclear shape redistributes the single particle states. The re- 
arrangement of the shells translates in modifications of the total deformation energy. 
These fluctuations, together with the influence of the pairing energy can be taken 
into account by means of the Strutinsky prescriptions [18, 19, 20]. In the so called 
macroscopic-microscopic model, oscillatory shell and pairing corrections were added 
to the smoothly varying liquid drop energy surface. The macroscopic-microscopic 
energy surface exhibits minims at hyperdeformations, especially where the ratio be- 
tween the length of the nucleus and its central diameter is about 2, reflecting a 
stability of shape isomers and the occurrence of deformed closed shells {21, 22]. 


The experimental mass distributions of fission fragments [23] revealed that the 
nuclear structure plays a primordial role in the fission process. It was suggested that 
the asymmetric mass division of fission products is due to a possible preference for 
N=82 magic configuration of the fragments [24]. A first qualitative explanation of 
the mass-asymmetry of the fission products was given in Ref. [25]. By searching 
the barrier saddle points of transuranium elements within the Srutinsky prescrip- 
tions, it was evidenced that the shapes characterizing the first barrier are reflection 
symmetric, while those corresponding to second barrier are mass-asymmetric. Reg- 
ularities in fission mass distribution were exploited in empirical models by relating 
them to simple mathematical functions [26]. Assuming a statistical equilibrium and 
involving only three parameters, that is a distance between the complementary frag- 
ments, an intrinsic excitation energy, and a collective temperature, correlated with 
the potential energies of the fragments, the trends of the of the isotopic yields in the 
mass and charge distributions, together with their collective kinetic energies were 
reproduced [27]. Later on, the mass and energy distributions of the fragments were 
also systematized as function of the primary fragment deformations in three main 
groups: standard I connected to the spherical shell N=82, the standard II connected 
with N=88 and the superlong channels for deformed fragments [28]. By extracting 
the weights of these fission modes from experimental data, it is possible to improve 
the evaluations of other quantities that depend on the mass distributions [29], as 
the neutron multiplicities described by sawtooth distributions. 


The mass and charge distributions of the fragments depends on the dynamics of 
the nuclear system. The fission fragments cannot be formed in their ground states 
due to energetic restrictions at the scission point, even in the case of adiabatic pro- 
cesses [30, 31]. According to the Bohr-Wheeler model, the fission process can be 
understand as being realized in two main steps. The first one treats the evolution 
of the nuclear system up to the top of the outer barrier. A slow adiabatic motion 
is assumed and the nuclear system behaves as a compound nucleus in thermal equi- 
librium. The second step consists in a rapid descent from the top of the barrier up 
to scission, the energy release being converted in internal excitation or translated 
in collective kinetic energy. The modality in which the energy is shared between 
the internal degrees of freedom and the macroscopic ones depends on the tensor of 
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inertia and of the viscosity. As discussed in Ref. [32], the dynamics of the process 
should explains the broadening of the mass and charge distributions as function of 
the damping. Two classes of theories are envisaged. The first one provides essen- 
tially deterministic solutions by means of the equations of motion in the collective 
degrees of freedom, the damping being small [33]. The second one considers a strong 
damping and use statistical arguments. In a dynamical theory of the first class, the 
boundary conditions at the saddle point determines uniquely the scission configura- 
tion, as obtained in Ref. [83] by mean of the Hamilton’s equations of motion. If an 
excitation is considered at the top of the barrier, it follows that a statistical distribu- 
tion of initial conditions can be envisaged, and the mass distribution is obtained by 
means of the Liouville theorem. In this formalism, an initial phase-space distribu- 
tion is connected with a final one. In the statistical theory, an equilibrium between 
all degrees of freedom is assumed during the descent of the barrier, that is a large 
dumping hypothesis. A density of states for the nascent fragments at scission can 
be calculated, being proportional to the mass distribution function. As mentioned 
Ref. [32], the two classes of theories become equivalent if an additional dumping can 
be simulated in the deterministic theories by taking into account mass asymmetry 
vibrations superimposed on the collective radial motion. These additional vibra- 
tions may be obtained by solving the time dependent Schrodinger equation in the 
mass-asymmetry coordinate as given by the fragmentation theory [34]. Although, 
it could explain the broadening of the mass distributions, the statistical model was 
not able to provide information about the dissipation until recently. By exploiting 
the strong dissipative coupling assumption and by considering the competition be- 
tween one-body dissipation forces and the cohesive ones, the fragment distributions 
and the most probable fragment kinetic energies were reproduced [35]. A precise 
measurement of isotopic distributions of fission products as function of initial exci- 
tation energy, by providing information about the kinetic energies, and the angular 
distributions will contribute to a better understanding of the dynamics of the pro- 
cess. Surprisingly good predictions of mass and kinetic energies distributions were 
obtained within the Langevin equations of motion [36], recently. 


The sawtooth shape of the neutron multiplicity [37, 38, 39] was also explained 
qualitatively by the asymmetric distribution of fission fragments and the structure 
of the deformed issued fragments [40, 41]. The neutron multiplicities were used to 
evaluate the excitation energies as function of the fragment mass [42], evidencing 
a higher excitation energy for the symmetric fragmentation. This behavior was 
considered as the evidence of the existence of at least two fission barriers [42, 38]: 
one for the symmetric fission that leads to excited fragments, and another for the 
asymmetric fission with moderately excited fragments and high collective kinetic 
energies. Later on, it was suggested that the neutron emission curves are mainly 
related to the structure of the deformed final fragments [41] and the excitation due to 
the dynamics [43]. The neutrons evaporated from fully accelerated fission fragments 
can be evaluated within statistical prescriptions. A review of the models used for 
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Figure 1: (a) Double fission barrier obtained by means of three smoothed joined 
parabola plotted as function of a dimensionless elongation variable named ¢«. The 
parameters of the fission barrier are taken from the work of Ref. [88] and corresponds 
the compound nucleus 7°U. The heights and the minimum of the barrier V are 5.7, 
2.2, and 5.7 MeV, with the corresponding stiffnesses hw 1, 1 and 0.6 MeV. (b) 
Dependence between the curve that gibes the logarithmic of the penetrability P and 
the excitation energy E* of the compound nucleus. The Schrodinger equation was 
solved within the exact method by extending the formalism of Ref. [75]. The 6- 
resonances in the penetrabilities are connected to the eigenstates constructed in the 
second potential well displayed in panel (a). The compound resonant states created 
in the first well are ‘filtered’ by the resonances of the second well. 
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such purposes can be found in Ref. [44]. A good agreement with experimental data 
was obtained in the framework of the Los Alamos model [45] which is based on 
the knowledge of experimental energy releases. Subsequent developments improved 
input parameters by means of multiparametric matrices as function of fragment 
masses, charges and total kinetic energies [46]. 


A delicate balance between of cohesive nuclear forces on one side, and Coulomb 
and rotational forces of disruptive character, on the other side, together with fluc- 
tuating shell effects and pairing effects gives rise to the deformation energy of the 
nucleus, that manifests as a surface in the configuration space spanned by the col- 
lective variables. For very large angular momentum, the fission barrier may even 
disappear [47]. Due to the rotational degree of freedom, an anisotropy in the angu- 
lar distribution of fission fragments is produced, being strongly related on the spin 
projection K of the transition states on the fission axis. This distribution is a direct 
consequence of the conservation of the angular momentum, as evidenced in Ref. 
[48]. For lighter nuclei, it was observed experimentally that the angular distribu- 
tions exhibit abrupt variations at energies close to the barrier top. These variations 
were interpreted as a manifestation of the role played by each fission channel. For 
heavier nuclei, the variation in the angular distribution is attenuated at threshold 
energies, the manifestation of the structure being attenuated. It was argued [49] 
that this behavior reveals the presence of an intermediate equilibrium state due to 
the existence of the isomer minimum. Accordingly, in the case of heavier nuclei, the 
transition states do not modify so drastically the angular distributions, their external 
barrier heights being lowered. The presence of an intermediate equilibrium allows 
the treatments in the frameworks of statistical approaches. A statistical model able 
to give information about the fragment angular distributions was deduced in Refs. 
(50, 49, 51, 52]. By using statistical approaches related to the scission configurations, 
an effective momentum of inertia is assumed constructed as a combination between 
the parallel and the perpendicular rigid body components. The model provide infor- 
mation about the effective momentum of inertia from the measurement of the mean 
square value of the fragments angular momenta projections Ko on the fission axis 
at the saddle configuration. Therefore some hints concerning the hyperdeformations 
are guessed. Based on this model, investigations based on experimental data [53] for 
transuranium elements indicated that the effective momentum of inertia is compat- 
ible to the outer barrier deformation at low excitation energies. For high excitation 
energies, the effective moment of inertia becomes compatible to calculations realized 
to smaller deformations, that are in accordance to saddle configurations calculated 
within the liquid drop model. These results show that the shell effects vanish with 
excitation energies [54] or the order of 50 MeV. For lighter elements, the deforma- 
tions calculated in this way by using the anisotropy data are in agreement with the 
liquid drop model expectations. The modifications of the angular distributions for 
different mass regions were in accordance with theoretical results. For example, a 
shift of the saddle point position to smaller deformations was predicted with the 
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two center shell model [55] from the mass region A=200 to A >230. In this model, 
the properties of the fragments manifest in the region of the outer barrier. For mass 
A 230, the fragments shell effects are important, as a influence of the strong shell 
effects of the doubly magic nucleus !°?Sn. It was also suggested that the variation of 
Ko at low excitation energy is due to persistence of superfluid effects [56, 57]. It is 
important to underline that an accurate measurement of the angular distributions 
can offer valuable information about the nuclear structure at hyperdeformation to 
test the validity of actual models. 


From the previous considerations, it follows that the distributions of the to- 
tal angular momenta of the primary fission fragments carry information about the 
whole dynamics and the fission mechanism. The initial angular momentum of the 
compound nucleus is redistributed not only on the angular momentum of the relative 
collective motion, but also between intrinsic fragments spins [58]. It was experimen- 
tally found that the average spin of the primary fragments in low energy fission is 
about 7 h [59]. However, the fragments spins have rapid fluctuations that depend 
essentially on the fragment mass and the mass asymmetry of the reaction. Unfor- 
tunately, there is currently no consistent theory to predicts the spin distributions 
in the fragments. A better experimental description will definitely contributes to a 
development in the theory. As already mentioned, the actual understanding of the 
angular distributions of fission fragments is based on statistical arguments, by taking 
into account an effective moment of inertia [50]. This approach was developed by 
including the dynamics of the fission [60], the latter achievements being described 
in Ref. [52]. 


Now, it is established that the fission barrier in the actinide region is charac- 
terized by a shape with two or even three maximums [61]. This barrier offers an 
explanation for the experimental findings, stating with fission isomers, structures in 
the fission cross section, or angular distributions. Such a barrier is displayed in Fig. 
1, approximated with parabolas, as done in general in evaluations. In the region 
of the threshold excitation energies, that is, close to the maximum of the barrier, 
the fission cross section reveals a structure formed by a large number of interme- 
diate resonances. These structures are caused by many discrete excited states that 
exist in the second well of the potential, called class II states. If the fission cross 
section is calculated as function of the excitation energy, the class II states give 
rise to resonances in the cross section. To infer the information about the fission 
shape, it is possible to measure the isomer excitation energy. It is now possible to 
identify the shape isomers, being obtained by a large variety of nuclear reactions if 
the energy induced in the system is sufficient to surpass the inner barrier. In this 
way, some patterns due to the process determined by a neutron inelastic scattering 
to the isomer, involving also a prescission neutron evaporation can be identified in 
the fission cross section [62]. Also, the rotational bands that are constructed in 
the second minimum offer information about the hyperdeformed nuclear states. It 
was found, a verified theoretically within the cranking model, that the rotational 
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parameter h?/(2J) is of the order of 3.5 keV, for isomer states. Accordingly, the 
sub-threshold fission exhibits cluster of resonances, due to the fact that the fission 
strength is modulated by the class II states. ; From experimenatl data, the coupling 
matrix elements between the class II states and the class I states or about the fis- 
sion of class II states can be obtained by an analyze of the detailed structure of 
each cluster. The averaged values of the coupling matrix elements offers directly 
the order of magnitude of the penetrabilities of the first and second barrier for a 
given excitation. Also, resonance doublets are identified [63]. They are caused by 
a degeneracy of the unperturbed class I and class II states. By using a polarized 
neutron beam, it is possible to assign unambiguously the spin to the detected res- 
onances at low energy [64]. Unfortunately, theoretical predictions of such resonant 
peaks cannot be realized realistically in terms of actual models, their investigation 
being concentrated to the activity of analyze. 


The scientific community was surprised by many features experimentally de- 
tected related to the nuclear fission process. These data developed the theoretical 
treatments, and step by step, the theory has become increasingly predictive. There- 
fore, by developing the shell models an island of stability for superheavy nuclei 
was predicted [65] or and by understanding the fission mechanisms the spontaneous 
heavy ion emission [66] is assessed from model calculations. In the theory evolution, 
many new experimental data benchmark the development. The fission still offer new 
behavior. For a recent example, a new type of asymmetric fission was experimen- 
tally observed in the neutron rich '8°Hg nucleus [67]. The theory expected that a 
symmetric distribution of fission fragments should be produced, with an emphasized 
yield related to the two semi-magic 9°Zr products. The explanation that followed in 
the framework of the macroscopic-microscopic approach improved the actual knowl- 
edge by revealing a so-called local minimum in the potential energy surface [68]. 
The recent developments of this subject are reviewed in Ref. [69]. The bimodal 
fission phenomena [70] is another recent discovery, in which two components of the 
kinetic energy spectrum of fission fragments are identified. Moreover, the inver- 
sion of the odd-even effect of the low energy fission fragment distributions observed 
experimentally is difficult to be reproduced theoretically. The odd-odd yields are 
larger for excitation energies of the fragments lower that 4-5 MeV [71, 72], an effect 
that can be characterized as strange and that cannot be explained statistically. It 
was understand by mean of a new dynamical pair breaking effect introduced in the 
microscopic equations of motion [73]. 


Despite the improved the theoretical description realized along the time, the fis- 
sion cross section is still difficult to be reproduced by using microscopic theories, 
and it is evaluated within phenomenological approaches. For instance, the maxi- 
mums of the double fission barrier are deduced in an empiric way by using some 
simple parametrizations for the shape of the barrier and a predefined nuclear level 
density. As stated many years ago, the cross section is simply proportional with the 
number of states obtained semi-empirically in the transient point or saddle point 
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configuration of the barrier following the Bohr-Wheeler hypothesis. The fission 
cross sections is evaluated always in accordance with some input parameters that 
rely on experimental observable. For example, at least 6 parameters are needed to 
reproduce approximately the shape of the fundamental fission barrier. These are, 3 
values for the extremes of the potential barrier plus 3 values of the corresponding 
stiffness parameters of each region of the barrier. In the frame of this approach, it 
was estimated in Ref. [74] that the fission barrier heights can be extracted with an 
uncertainty of 0.5-1.0 MeV. Moreover, if the resonant structure of the fission cross 
section should be reproduced, then the situation is more ambiguous because many 
transition states are introduced empirically, by hand, to fit the experimental data. 
So, it can be concluded that this procedure does not include the resonances due to 
the beta vibrational and to the rotational couplings adequately, in order to under- 
stand microscopically the fission process. Therefore, fission of actinides as well as 
sub-actinides is not yet understood sufficiently for increasing energies above a few 
MeV. The understanding of the experimental findings still require not only accurate 
experimental inputs but a sustained effort for their theoretical explanations. At this 
point, the phenomenological approach that is used in evaluations for the fission cross 
sections analysis is still limited by many ambiguities and shortcuts. Some of them 
are mentioned in the following. 

In evaluations, the double or triple fission barriers are represented with several 
smoother joined parabolas [75]. The deviation of the realistic barrier shape from the 
harmonic one is a cause for an alteration of the values of the penetrabilities. In the 
approximation used, only at energies close to the fission barrier maximum, confident 
values of the penetrabilities can be calculated. 


An average of the fission cross section is simulated by using a set of values for the 
heights of the barriers and for the nuclear level density. Unfortunately, it is known 
that two different sets of values for these quantities can give the same value of the 
cross section in special circumstances. 


The nuclear densities depend on the deformation of the parent nucleus and on 
it fission trajectory [76]. For the hyperdeformations compatible the region of the 
outer barrier, the level density is predicted to increase to very large values. One of 
the reason is the fact that the momentum of inertia is enhanced at hyperdeforma- 
tions. Therefore the number of transient states is increased. In the evaluations, the 
variation of the level densities during the fission is treated mainly with phenomeno- 
logical corrections, some of them in the aim of fitting the data. One of them, the 
Porter-Thomas distributions [77] are based on statistical analysis of the fluctuations 
of the reduced widths. Collective vibrational and rotational enhancements of nu- 
clear level densities close to the top of fission barriers in principle can be evaluated 
microscopically but are also considered in a phenomenological way. 

Unfortunately the phenomenological approach is considered as valid for predic- 
tions. For example, the heights of the fission barriers obtained from fission data are 
gathered in systematic and used afterward to predict the cross section of unmeasured 
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nuclei. In this philosophy, the nuclear structure, that may change dramatically from 
one nucleus to another, plays a minor role. But, it was seen experimentally that one 
neutron can change the overall mass distribution from asymmetric to symmetric. A 
such result is reported in the case of 7°%?°$Fm [78], the existence of at least two 
competing fission trajectories is postulated. 


If a triple fission barrier is postulated, the evaluations show that a dine struc- 
ture is produced in the fission cross section resonances at the threshold energies, 
phenomenon generically known as the Th anomaly [79, 80]. A shallow ternary min- 
imum of about 1 MeV was predicted within the macroscopic-microscopic approach 
at a large amplitude asymmetric shape [81]. As given by the theory, the mass- 
asymmetry parameter related the third minimum is much more evident than that of 
the second minimum or that that of the scission point. If the least action principle 
is used to characterize a fission path it is difficult to fit this behavior. The fission 
trajectory has to go through different regions of the configuration space to feel the 
third well. But, a very large increase of the inertia is produced in the inflection point 
of the trajectory decreasing tremendously the penetrability of the process. Another 
way to attack the Th anomaly is with dynamical single particle effects [82]. 


The variation of the nuclear inertia is neglected. It could be incorporated by 
fitting the widths of the parabolic potential barriers. The stiffness parameters of the 
parabolas influence the widths of the evaluated cross section resonances. Variation of 
the nuclear inertia can produce supra-barrier resonances [83]. Moreover, the inertia 
varies also with the excitation energy [84]. In a realistic treatment, not only the 
widths of the resonances can be degraded, but their positions are altered. 


A nuclear density is a statistical parameter can be calculated only for large values 
of the excitation energies, where the number of nuclear levels per energy is large. For 
excitation smaller than several MeV, the excitation spectrum is essentially discrete 
and can be reproduced with the help of a deformed shell model. But, in the actual 
calculations even the low energy transition states are evaluated empirically in terms 
of densities to participate in the fission cross section. However, they are calculations 
in which some transition states are introduced empirically in the calculations with 
the aim to reproduce the positions of the experimental fission cross section reso- 
nances. These transition states are treated simply as a sequence of levels ordered by 
their spin that cannot usually intersect during the deformation of the parent nucleus, 
indexed with different values of the excitation energies in the extreme of the barrier. 
The number of fitting parameters increases, because for each added excitation the 
number of input parameters increases with 3, that is, the values of the excitation 
at the extreme of the barrier. These states are kept in the same ordering, and the 
existence of the avoided levels crossing region is completely neglected. Especially 
in the region of the outer barrier, the calculations showed that the density of sin- 
gle particle levels increases, and the probability to produce avoided levels crossings 
regions is very high. This effect is also reflected by the values of the shell effects, 
being decreased. It is also known that when a prolate nucleus deforms in its way 
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to scission, the orbits of the high spin states are preferentially located in a plane 
perpendicular on the elongation [63]. These orbits are forced to move into a smaller 
restricted value of the diameter of the nucleus when the nucleus deforms. The single 
particle energies of the levels with high intrinsic spin should increase. Accordingly, 
a lot of crossings between levels should be produced during the dynamical evolution, 
the levels being rearranged. A microscopic calculation should be realized prior to 
the evaluation in order to identify such level crossings. 


The population of the transition states is considered essentially unity. Such 
assumption is valid only in the case of an adiabatic system, by neglecting residual 
interactions, as the pairing for example. But, the probability of occupation of a 
transition state should depend on the dynamics of the process [85]. Due to the 
rearrangement of the single particle levels, the probabilities of occupation of these 
states are strongly influenced also by the inherent radial and angular couplings 
between single particle levels. The probability of occupation increases or decreases 
as function of the residual interactions, being not compatible to that of an adiabatic 
system. Accordingly, the population of the transition states should be obtained by 
solving the microscopic equations of motion. 


Due to an incomplete knowledge of the fission mechanism, sometimes a complex 
fission barrier is used to figure an effective damping, to take into account in a 
phenomenological way additional de-excitation processes. It is expected that by 
using an imaginary component of the barrier it is possible to simulate the absorption 
of the fission degree of freedom into more complicated states caused by additional 
interactions. As possible processes one can note the y-deexcitation in the second 
minimum, or the emission of one pre-scission neutron. However, the description of 
the coupling between the collective motion and the internal degrees of freedom to 
infer the dissipation is not possible within this model. By introducing an imaginary 
component of the barrier in the second well, one obtains two effects: the transmission 
is reduced and broader resonances are simulated. Unfortunately, the dissipation is 
not evaluated. Because in the optical model, the fission fragments resort without 
internal excitation. 


In neutron induced fission, a large excitation energy is brought to the compound 
nucleus. This energy is evaluated in a statistical way, by considering that it is 
equiprobably shared between intrinsic and external degrees of freedom, as usually 
considered in the ergotic theory. That is, the collective kinetic energy of the nuclear 
system can have any value comprised between zero and the maximal amount of the 
excitation. As mentioned, only one set of parameters characterizes the fission barrier 
in wide interval of intrinsic excitation energies of the compound nucleus. The cross 
section fission at a given excitation energy is obtained as a integral over the excitation 
energy, by spanning with the same probability all the possible values. In these 
circumstances, for a correct evaluation, the parameters of the fission barrier should 
be modified when the intrinsic excitation energy increases, due to the disappearance 
of the shell effects [86]. Always, we are left with a distribution of fission barriers 
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that depend on the excitation energy. The parameters of the fission barrier should 
vary in accordance with the intrinsic excitation energy. The actual evaluations 
cannot include a such phenomenon in calculations, because they need a theoretical 
calculation of the shell effects at each deformation. The presumed deformation of 
the parametrized fission barrier is not consistent to a real deformation associated to 
a degree of freedom obtained from a theoretical model. 


The models used in evaluations require the knowledge of the fragments excita- 
tion energies. But these information are usually extracted from the experimental 
data, that is from the quantities that should be theoretically evaluated. The fission 
cross section data are very complex and carry information about all the interdepen- 
dent processes, about the nuclear structure, an more important, about effects due to 
the dynamics of the fragmentation. A very important observable is the dissipation 
energy obtained at the nuclear rupture. From the theoretical point of view, infor- 
mation about the dissipation can be obtained only from a dynamical treatment of 
the fission process. But it is difficult to unfold the dynamical effect from the puzzle 
created by different concurrent mechanisms. 


Concerning the evaluations for fission cross sections purposes, it should be reiter- 
ated that they depends always on some input parameters extracted from experimen- 
tal data. So, a vicious circle is created, the evaluated data are based on correlations 
derived from experimental findings and will not be able to exceed the accuracy with 
which the findings were obtained [87]. Consequently, the evaluations used in the 
design of nuclear reactors and other applications demand a high accuracy of exper- 
imental data in order to provide confident results. But, in the same time, these 
evaluations cannot be used confidently to make predictions. Some of the previous 
ambiguities are know from long time [88, 89, 90], and summarized recently in Ref. 
[91], but no satisfactory solutions were found. Therefore, if a such analyze is applied 
to a specific isotope, it provides information about some observable that are only 
appropriate to the investigated nucleus. The predictions concerning the unmeasured 
nuclei are questionable due to the lack of consistency. These peculiar aspects can be 
only solved by investigating microscopically the fission process aiming to provide new 
information about its basic mechanism [92]. The microscopic approaches used in the 
fission investigations are subject to many constraints in order to be credible. They 
should be able to reproduce a large set of experimental data: fission barriers, fission 
mass distributions, total kinetic energies of fission-fragments, spontaneous lifetimes. 
It is a huge task. As we will see, theoretical models have yet to be improved to 
provide data with consistencies comparable to those obtained in fission cross section 
experiments. For the moment, the best way to solve the enumerated problems is 
to improve our actual experimental knowledge. The experimental data carry in- 
formation about all the interdependent mechanism that follow from the structure 
and the dynamics of fission. Therefore, they should be used as benchmarks for the 
theoretical investigations. 


The prediction degree of the evaluations is strongly related on the accuracy of 
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the experimental data as well as on the way in which the theoretical achievements 
are incorporated. Progresses were realized in the statistical modeling of the fission 
observable corroborating as much as possible information obtained from theoretical 
concepts [93, 94] or general properties of the nuclear matter. The model predicts 
the Thomas-Fermi fission barriers in accordance to the topographic theorem [95] 
by including final shell effects, the individualization of the fission fragments being 
dictated by the two center shell model behavior [55]. A large spectra of fission 
observable can be evaluated. Recently, the model was improved by postulating a 
sorting of the excitation energy at scission [96] modulated by the pairing energy gaps 
of the fragments. By taking into account the differences in the nuclear densities of 
the fission partners, final excitations of the order of 3-8 MeV are obtained. 


The basic ingredient used in the theoretical study of the large amplitude motion 
processes is the deformation energy of the nucleus. Apart the phenomenological 
approaches, macroscopic-microscopic models or Hartree-Fock ones are used to ap- 
proach the fission process [97, 98]. Within the macroscopic-microscopic method, the 
measured ground-state energies were reproduced with deviations of 1 MeV [99, 100]. 
Instead, the values of the fission barrier must be obtained with almost the same 
uncertainties [101, 102]. The nuclear shapes around scission needs many degrees 
of freedom to be characterized conveniently [103]. In order to rend the problem 
tractable, we limit ourselves to a small number of generalized coordinates, that is 
the elongation, the necking, the deformation of the fragments, the mass-asymmetry, 
and axial y-asymmetry. These drawbacks leads to an insufficient description of 
hyperdeformations, or scission configurations. Moreover, when a parent nucleus un- 
dergoes fission, in the scission point, the problem of individualization of the fission 
fragments is still unsolved. In the scission configuration a unique entangled distri- 
bution of microscopic states is calculated. The numbers of proton and of neutron in 
the two fission fragments are obtained as averages. 


In the macroscopic-microscopic approach [104, 105], a parametrization of the in- 
dependent particle potential should be supplied. There are no arguments to consider 
this parametrization valid for nuclei far from the line of stability or for hyperdefor- 
mations. A more rigorous alternative is to rely on self-consistent potentials given 
by effective interactions between nucleons which are not related to deformations 
(106, 107, 108, 109, 110, 111]. In the Hartree-Fock approximation, a smooth self- 
consistent mean field is obtained because it includes the convolution of many instan- 
taneous densities. If this mean field is varying in time, due to the modifications of 
some constraints, each single particle wave function is moving independently in the 
mean field. The Pauli principle is fulfilled. The single particle collision, that must 
be the cause of the mean field, is incorporated in the equations of motion only in the 
measure in which they contribute to the mean field potential. Usually, a state depen- 
dent pairing model is also used [112]. Unfortunately, the pairing mean field is still 
considered as function of some parameters as in the density dependent delta interac- 
tion approach [113, 114, 115]. The dynamics of the fission process depends strongly 
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on the parametrization used for the pairing mean field [112, 116]. As mentioned 
above for the macroscopic-microscopic models, it is questionable if the parametriza- 
tion is still valid at hyperdeformations. The phenomenological assessed barriers 
represent a test for the effective interactions. The fission barriers calculated with 
Hartree-Fock models, using different effective interactions [117, 118, 119, 120, 121] 
show that the mean field models have become competitive with the macroscopic- 
microscopic ones. 

For large scale collective motion, in general, two approaches are used. On one 
hand, the generator coordinate method assumes that the internal structure of the 
evolving systems is equilibrated at each step of the collective motion [122]. On the 
other hand, the exchange between the collective and intrinsic degrees of freedom 
is neglected, that is, adiabaticity is assumed [123, 124, 125]. In the low energy 
domain, when the velocity of the nucleons in the mean field are much greater than 
the velocities for the modifications of the generalized coordinates, the dynamics 
problem can be solved using equations of motion or quantum statistics. Up to now, a 
clear response was not given concerning the modality in which the dissipation energy 
intervenes in a quantum assembly. In the quantum statistic formalism, it is possible 
to introduce concepts as friction or viscosity [126], in analogy with the classical 
mechanics. Concerning the microscopic equations of motion, the problem remains 
ambiguous. The first step to be overcome is to write down of some time dependent 
equations for the intrinsic variables. To do that, the variational principle can be used 
or the Heinsenberg form of the equations of motion [127, 128]. In the last kind of 
treatment, a complete solution that describes the dynamics of the system is obtained. 
This solution reflects better the response of the nuclear system to the changes of 
the single particle mean potential. But, there is no clear answer to how to get the 
dissipation from these results [129]. Dissipation is an irreversible flow of energy 
and angular momenta from collective degrees of freedom into intrinsic ones. In this 
sense, both macroscopic-microscopic and Hartree-Fock dynamical models should be 
improved to be able to reproduce the experimental findings. It should be noticed 
that uncertainties of only one MeV in the fission barrier heights can be translated 
in deviations of many order of magnitudes in the calculation of spontaneous fission 
the half-lives. 

The tendency today is to use the Hartree-Fock models. Because, it is assumed 
that starting from fundamental interactions it will be possible to explain as a whole 
most of the mysteries surrounding the fission process. The models used in present 
are able to reproduce only some facets of the richness of the process behavior, and 
cannot be extended to predict other encountered phenomena. 


3 Dynamics 


Nuclear reactions were intensively studied in the last thirty years, one of the most 
intriguing aspects being that the nuclear dynamics is lead in a regime in which the 
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statistical physics has an important role. Therefore, concepts as friction or thermal- 
ization were introduced in this field. Actually, advanced investigations are ongoing 
to test the basis of these concepts when a small and dense number of fermions are 
involved, as in the atomic nucleus. The majority of the applications cannot be 
treated within simple thermodynamic models and require sophisticated concepts of 
non-equilibrium statistical mechanics, especially to take into consideration of the 
strong fluctuations. In this sense, as an example we enumerate the fission of ultra 
hot nuclei and the collisions of heavy ions that lead to very excited nuclei together 
with multifragmentation. The development of some improved approaches is also 
motivated by the recent availability of some ion beam facilities ranging almost all 
the energy domains, from Fermi energies up to relativistic ones. A great part of the 
theories dealing within these aspects borrowed the idea of non-equilibrium statistic 
physics, especially the kinetic theory and the stochastic method. These approaches 
were extended and adapted to the particularities of the nuclear dynamics to treat 
strong correlated fermions. These theories represent a critical test for the kinetic 
equations in a strong dissipative regime that implies fluctuations. 


In the low energy domain, when the velocity of the nucleons in the mean field 
are much greater than the velocities for the modifications of the generalized coor- 
dinates, the quantum statistic intervene. One behavior observed experimentally is 
that the fission fragments are usually strongly excited. That is, the process cannot 
be adiabatic. A flow of energy and angular momentum is produced during the decay, 
from collective degree of freedom to intrinsic ones. Up to now, a clear response was 
not given to the modality in which the dissipation energy intervenes in a quantum 
assembly. The first step to be overcome is to write down of sometime dependent 
equations for the intrinsic variables. To do that, the variational principle can be 
used or the Heinsenberg form of the equations of motion [127, 128]. In the last 
kind of treatment, a complete solution that describes the dynamics of the system 
is obtained. This solution reflects better the response of the nuclear system to the 
changes of the single particle mean potential. The same solution is obtained within 
both methods only if the antisymmetric matrix elements of the time derivative op- 
erator are neglected. Using this last approximation, it is believed that the major 
part of the collective energy associated with the coherent motion of the nucleons is 
eliminated. This part,which is in general neglected, must define the effective mass 
tensor of the collective motion together with the its dependence of the excitation 
energy. 


The dynamics of the nuclear fission was investigated by solving the microscopic 
equations of motion in Refs. [132, 133, 82, 134, 73, 135, 136, 137, 139, 84, 140] The 
problem of diabaticity was explored by using a monopole pairing force, sufficiently 
weak to no strongly redistribute the nucleons so that the mean field does not change 
significantly. It should be noted that the study available in reference [114] assessed 
that the seniority state independent pairing force begins to be insufficiently accurate 
only for nuclei close to the stability lines that characterize the spontaneous emissions 
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of protons and neutrons in the nuclear map. Close to the stability valley, the BCS 
approximation used in this project is performing well enough as it can be seen in 
the global fit of mass realized in reference [130] or in the comparison concerning the 
density dependent delta interactions in the seniority force in the case of occupation 
probabilities of single particle levels [131]. 


In a normal way, the nuclear system is parametrized with several generalized 
coordinates that parametrize the shape. The independent degrees of freedom are 
associated to these coordinates. By mean of these generalized coordinates the be- 
havior of all other intrinsic coordinates are determined. The basic ingredient for 
such an analyze is a nuclear shape parametrization that depends of the degrees of 
freedom. The generalized variables associated to the degrees of freedom vary in 
time leading to a compound system and to the split of the nuclear system in two 
separated fragments. A microscopic potential can be constructed to be consistent 
within this nuclear shape parametrization. It is known that a nuclear shape can be 
used for binary processes if the following conditions are satisfied: (I) the three most 
important degrees of freedom [148], that is elongation, necking and mass asymmetry 
are taken into consideration; (ii) A sphere and two separated fragments are allowed 
configurations; (iii) The smoothness of the neck is a independent variable. These 
three conditions are satisfied by the parametrization used in the following: two el- 
lipsoids with different eccentricities and semi-axis smoothly joined with a median 
region given by an arc of circle rotated around the axis of symmetry. A good model 
that can be used in nuclear fission must have at least five degrees of freedom, that 
is the elongation, the necking, the mass asymmetry and the deformations of the 
two nascent fragments formed during fission. As our model, this parametrization 
can characterizes both the configuration associated to a single fragment and that of 
two separated ones. Also, for this nuclear shape parametrization we will use a two 
center Woods-Saxon shell model realized by the authors of this proposal. Using this 
model we are able to solve the mean field potential so that to have only one system 
of orthogonal wave functions for a two body system that approximate correctly the 
mean field. Because we have only one system of wave function it is possible to solve 
correctly the antisymmetrization problem of the total wave function. In the case of 
one center potentials, the problem of two bodies is solved using molecular methods, 
using two separated potentials, and the antisymmetry of the total wave function 
cannot be obtained exactly. Actually, we are the sole group able to use two center 
Woods-Saxon potentials. 


To determine the deformation energy of the system, we use the macroscopic- 
microscopic method, that is, the total energy is given by the sum of the liquid drop 
energy and the Strutinski shell and pairing corrections. The liquid drop energy will 
be computed in the frame of the Yukawa plus exponential model extended for bi- 
nary systems with different charge densities. The diffusivity electrostatic energy, the 
Coulomb energy, the surface energy, the volume energy are all included together. 
To compute the microscopic shell correction we use the single particle diagrams ob- 
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Figure 2: Schematic representation of an avoided levels crossing region. The Landau- 
Zener effect is produced when a such avoided levels crossing region is produced during 
the deformation of the nuclear system. 


tained within the two center Woods-Saxon shell model.The nuclear inertia will be 
computed in the frame of the cranking model or in the frame of the gaussian overlap 
approximation. Once we have a Hamiltonian that contains also the pairing interac- 
tions and we introduce the collective parameters through the Lagrange multipliers, 
the response of the nuclear system for shape changes is obtained with the effective 
cranking mass. A complex code for this model in the adiabatic hypothesis have been 
realized recently [135] for the two center Woods Saxon model. Therefore, the two 
main ingredients too calculate the action are at our disposal. 

The shape of the potential barrier can be obtained only if the trajectory of the 
nuclear system in the configuration space spanned by the five generalized coordinates 
associated to the degrees of freedom is known. This trajectory is obtained by the 
minimization of the action integral or by solving the Euler Lagrange equations of mo- 
tion. Up to now, the minimization methods were realized only for adiabatic systems. 
We will introduce in the frame of the method the calculation of the dissipated en- 
ergy within the time dependent pairing equations, that are similar to Hartree-Fock- 
Bogoliubov equations of motion. These equations are deduced from the variational 
method. The variation versus the time of the Bogoliubov wave function amplitudes 
in seniority 0, 1 or 2 states are obtained. The equations of motion for the mixing be- 
tween seniority 1 configurations will be deduced by incorporating the Landau-Zener 
dynamical effect and the Coriolis couplings, as generalized recently [141, 73]. 


3.1 The microscopic equations of motion for seniority-1 mixing 


Two microscopic effects participate to the mixing of the configurations during the 
large amplitude motion: the Landau-Zener promotion mechanism [142, 143, 144] 
and the Coriolis coupling [145]. Two levels with the same good quantum numbers 
associated to some symmetries of the nuclear system cannot intersect during the 
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deformation and exhibit avoided levels crossings regions. In Fig. 2, a such avoided 
levels crossing region is schematically represented. In such a region, the adiabatic 
levels €, and €g exchange their characteristics. If a nucleon is located initially on the 
adiabatic level €;, and the passage of the avoided levels crossing region is made slowly, 
the nucleon will remain on the same adiabatic state €; with a great probability. If 
the passage is made rapidly, with a large deformation velocity, the the nucleon will 
skip with a great probability on the superior adiabatic level €2. In this case, the 
nucleon follows the diabatic state e;. The interaction energy between the diabatic 
states e; and eg is given by the energy difference at the distance of closed approach 
between the two single particle levels, that is 2h19. 

In this section, as an example, the equations of motion for the seniority-1 nuclear 
systems are deduced. These differential equations takes into account the microscopic 
effects that cause the configuration mixing. In the framework of the variational 
principle, one has to minimize the functional 


6L=6(prm|H + Hp—-ihf + H! —A\N2N, — N\Np| 


erm) (1.1) 


where H is the many-body Hamiltonian which includes the pairing | interaction, 
Hp is the centrifugal energy, H’ is the term which allows the Landau-Zener effect 
[141], while Nj, and No are particles numbers operators acting on the two fragments 
that are formed in the scission region. The condition introduced in term of particle 
number operators determines a dynamical projection of the numbers of particles 
on the two final fragments [139]. Here, \ is a Lagrange multiplier. In the energy 
functional (1.1), the trial many body wave function is given as a superposition of 
seniority-1 Bogoliubov wave functions, allowing rotations. The trial wave function 
is defined as follows: 


lem) => combi vom I] (wor 1m) als Var, 1(0.m) Ay 4%) |0) (1.2) 
Q,m (2 ,1)A(Q,m) 


where CQ,mm are amplitudes. The square of these amplitudes supply the probability 
to obtain a given configuration. The creation and annihilation operators for the 
single particle state (2,1) are ag | and ag, respectively. The BCS occupation and 
vacancy amplitudes are vg, and ug, respectively. With 


Q = —|Q| (1.3) 


are denoted the intrinsic spin projections on the axis of symmetry. With 


GT a 0 ac 
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one denoted a product between the creation operator a8 mm of the blocked level and 
the rotation function Dig (w) . The relevant properties of the rotation function are 


— 
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being useful in the calculation of matrix elements. 

To deduce the equations of motion, first of all, the elements of the Hamiltonian 
should be obtained. For the single particle energies, after some algebra one obtains 
the next formula: 


< ap [[ - vay a7) | So ex(aj ap + az az) | af [[ - vay az ) > (1.8) 


I¢R k I¢R 
=2 a PLEL + Ex. 
It¢k 
where ¢; are single particle energies. The single particle densities py = |v|? and the 


pairing moment components «; = ut; are used in the final expressions. For the 
pairing rezidual interactions, one finds: 
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= ve — G isk pi 


where Ay = G)> ugvz is known as the pairing gap parameter. For the condition of 
normalization one obtains: 


(1.9) 
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One needs also to calculated the matrix elements of the expressions obtained 
with the angular momentum ladder operators. 


(yp lj+|y) = > in cr em (0 Jay Leen’ (wm) + Vion war ) j+ 
x lax, ie (re + vumyat af | 0) 
To effectuate this calculation, one identifies two cases. The first situation is for 


m’ = m, that is the matrix elements between the same seniority-1 configuration. 
That gives: 


(1.11) 
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All the previous terms should be zero. The second case is when the matrix element 
(Olam J+a;* |0) is different from zero. The next equations are obtained 
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The first term in the right hand side can be readily transformed as 
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By indexing with 2 and m the result is: 
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For n = m’, the second term in the right hand side gives: 
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By indexing with 2 and m, this term becomes: 
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The final results is: 
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That is, indexing with 2 and m, the matrix elements are 
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The matrix elements within the rotation functions are given by the next expression: 
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The rotational energy for the axial symmetric rotor is defined as: 
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Therefore, the matrix elements to be calculated are: 
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for Q = —I,..,/. It is also important to determine the next expressions: 
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(1.23) 
where the angular momenta j, = (j4 + j—)/2 and jy = (j4+ — j_)/2i are given in 
terms of ladder operators. That means also that when the term j? is evaluated, the 
terms of the type 7+ < j— are positive for ne In the same time, the terms of the 
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type jy X jx, j- X j- will be canceled with terms obtained from j2. The firs step is 
to calculate 
(Pm! |J+ + J-lPm) = Olam Ty gm! pe a Ui (my Wo) 
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Here, one has many possibilities: m’ = m where the integral is zero; m’ 4 m where 
the integral is zero with only one exception m’ = | and I’! = m In this last case, one 
obtains 


(Pm 3+ + 5-|m) 
= (Dlai(umay + UF gy) 4mam)(J+ + J-)am(Ui(m) + Vi(mn) a )10) (1.25) 
= Um (1) Un(m) Olar G+ + J—)amn|0) + VF, Pen) (lan he az |0) 


Therefore, it follows: 
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If the two previous expressions are summed together, one finds 
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[day lila) at lizlar) + (af li+laz) (art |j_-lajr)] 


(1.29) 


The completion relation is realized through the summation over the wave func- 
tions of all the intrinsic states states at, but is not directed by the wave functions 
b+ which includes the rotational functions D/,,(w). Finally, the relations read by 
indexing with the numbers 2, m as 
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(1.30) 
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The expected values of the energy functional is then: 
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(1.31) 
In order to minimize the functional, the expression (1.31) is derived with respect 
to the independent variables vj,,) and Uo) by taking into account the subsidiary 
condition uz + ve=1. Two relations follows, one being 


Yn | Gm FF 4 20ftyn (et — ren) 
Ym) Ym) Pl(m) 
[pn (~HLES) + (oy BBL) 
x Dagm “(my + 200m) YiGny| + *R8j om) f = 9 


and the other being its complex conjugate. This system can be solved by considering 
that the expression in the curly bracket is zero for each value of m. The time- 
dependent pairing equations associated to an unpaired nucleon in the state m emerge 
if no configuration mixings are allowed [146, 147]: 
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1.34 
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Now, within the same expected value of the functional, it is possible to obtain 
the equations for the mixing of configurations. These equations are similar to the 
equations needed to obtain the probability that an unpaired nucleon is located on 
a state m. For this purpose, the expression (1.31) must be derived with respect cj 
and ct,. Two equations follow, one being: 
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and the other being the complex conjugate. In the previous expression, the time 
derivative of Ug’ m(Qjm) intervenes. This derivative can be calculated by using the 
time-dependent pairing equations (1.32). The result is: 
(US) mn! (Om) 0M! mn! (2m) — BO) mm! (Q,m)YO!,m!(O,m)) 
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(1.35) 

Therefore, the time-dependent equations for configuration mixing is obtained: 
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If the amplitude of a seniority-1 configuration is known, and the nuclear system 
evolves in time, the previous equations give us the amplitudes of all other seniority- 
1 configurations that are populated during the evolution of the system. All these 
seniority-1 configurations are connected by “gates” opened by the Landau-Zener and 
the Coriolis interactions. Of course, another system of equations is available for the 
complex conjugates of the amplitudes.. 


In general, in nuclear physics it is usual to speak about probabilities instead of 
amplitudes. The equations can be recasted to gives the probabilities instead of the 
amplitudes. Multiplying the Eq. (1.36) with c,,, its complex conjugates with c*,, 
and subtracting, the next relation follows. 
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iFrom relations (1.36) and its complex conjugate, another relation can be de- 
duced 
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The derivatives / and 0* are still kept in the previous expression to make the equa- 
tions simpler. The Rel. (1.35) which is deduced from the time dependent pairing 
should be also recasted in terms of probabilities. The transformed relation follows: 
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This last relation reveals the fact that the equations for configuration mixing depends 
not only on the interactions in the avoided crossing region and on the Coriolis 
coupling, but also on the dynamical occupation probabilities that are obtained for 
the same collective velocity. It is interesting to note that the expressions inferring 
the time derivatives of the BCS occupation amplitudes in Eq. (1.38) have energy 
dimensions. 


The Rel (1.37) gives us the probabilities to find a seniority-1 configuration at a 
given deformation, for a given collective velocity. The second ecuation that can be 
deduced represent a moment component between different seniority-1 configurations: 
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After some rearrangements of terms, and defining the moment components of 
the probabilities $(Q,m)(Q1,m1) = CQ,mCO1,m, and transforming the Rels. (1.37) and 
(1.39), the Eqs. (1.40) and (1.41) are eventually obtained. These equations de- 
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pend on the single particle densities and the pairing moment components. So we 
obtain the dependence in time of the seniority-1 probabilities and of their mixing 
components. The equations follows: 
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These equations should be solved numerically by knowing the single particle states, 
the Landau-Zener and the Coriolis couplings. Simple examples with solutions of 
these generalized time dependent pairing equations are given in the following section. 


4 Results 


This section addresses several examples concerning the applications of the previous 
relations in fission dynamics. The first example concerns the fission of 2°?Th as 
investigated in Ref. [159]. By using the time-dependent pairing equations given by 
Rels. (1.33) and (1.34), it is possible to determine the fission times. This system 
of coupled equations are able to infer the excitation energy at scission as function 
of the internuclear velocity between the two nascent fragments. By comparing the 
theoretical values of the dissipated energies and the experimental ones, the best 
value of the internuclear velocity is obtained. By knowing the collective velocity 
and the shape of the barrier, the determination of the fission time is determined in 
an obvious way. 


But, in order to solve the system of the equations of motion, the single particle 
levels diagrams are required for the protons and neutrons configurations. That is, a 
disintegration path, the single particle diagrams, and some models to calculate the 
inertia, are required. First of all we need to know how the single particle levels are 
rearranged during fission, how they react to the change of the mean field. Usually, 
it is considered that the nuclear mean field is managed by some generalized coor- 
dinates, which correspond to the most important degrees of freedom of the system. 
The generalized coordinates vary in time leading to the scission. So, a nuclear shape 
parametrization is the main ingredient needed in our calculations. In the examples 
given in this work, the nuclear shape parametrization is axial symmetric, being a 
combination of two ellipsoids of different semi-axis smoothly joined by an interme- 
diate surface given by the rotation of an arc of circle around the axis of symmetry. 
This nuclear shape parametrization depends on five degrees of freedom: the elon- 
gation, the mass-asymmetry, the necking, the deformations of both fragments. The 
elongation is denoted R, being the distance between the centers of the spheroids. 
The mass asymmetry is denoted a;/ag, that is the ratio of the major semi axis of the 
two nascent fragments. The necking is represented by C3, the curvature of the arc 
of circle that defines the intermediate surface. The deformations of the fragments 
are b;/a; and b2/a2, that is the ratio between the semi-axis of the fragments 1 and 
2, respectively. The essential problem is to find a trajectory in the configuration 
space spanned by our five independent generalized coordinates that starts from the 
ground state and arrives at the scission, where the two complementary fragments are 
formed. For this purpose, one determines numerically a path in the configuration 
space within the algorithm proposed in Ref. [149], in which the dependencies of the 
parameters a;/a2, C3, b;/a, and b2/az are obtained as function of R. This path 


122 M. Mirea 


should minimize the penetrability given by the WKB formula 


Paesp{= | VV BIB Rar} (1.42) 


where V(R) is the deformation energy, B(R) is the inertia and Ep is the zero point 
kinetic energy. The macroscopic-microscopic model [150, 151] is used to calculate 
the deformation energy. In this model, a combination between a liquid drop energy 
and shell effect offers a good estimation o the total energy of the nucleus as function 
of the deformation parameters. To calculate the liquid drop part, several contribu- 
tions are taken into account: the Coulomb energy, the surface term which is given 
by the Yukawa-plus-exponential interaction [152], the Coulomb diffuseness correc- 
tion contribution, the volume term, and the Wigner one. The formulas for these 
contributions entering in the calculation of the liquid drop energy are presented in 
Ref. [153]. At scission, the model is generalized for binary systems with different 
charge densities [154]. The shell and pairing effects constituting the microscopic 
corrections are calculated by using the Strutinsky procedure [155]. By using this 
method, the rapid fluctuating effects due to the intrinsic structure can be inferred 
in the total energy of the nuclear system. The fissioning nuclear system acquires a 
potential that has a double humped shape. The intermediate well is responsible for 
the resonances that appear in the structure of the fission cross section, as displayed 
in Fig. 1 and explained in the text. The inertia along the least action trajectory is 


Oqi Og 
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where q; and q2 are two generalized coordinates, and By,q. is the element of the 
mass tensor versus g; and qo. This symmetric tensor describes the reaction of the 
nuclear system towards the external forces which modify the generalized coordinates 
q;. The tensor of inertia can be obtained from the variational principle as given in 
Refs. [84, 135, 140]. In Ref. [140], the behavior of the inertia for the nucleus 
investigated in this example is described. The main generalized coordinate in our 
model is the elongation (or the internuclear distance) R defined as the distance 
between the centers of the two nascent fragments. 

Once the fission path is obtained in the configuration space it is possible to 
calculated the single particle diagrams required to solve the equations of motion 
given by Rels. (1.33) and (1.34). The single particle levels were calculated within the 
Woods-Saxon two center shell model [141]. This model is able to describe the passage 
from one nucleus into two separated bodies in a very accurately way, being useful 
to characterize scission configurations [156, 157, 158]. Even more, this precision 
also holds true for very large mass-asymmetries, as obtained in the case of alpha 
or cluster decay [160, 161, 162, 163]. In our model, the ground state of ??4Th is 
deformed, the fundamental elongation being R ~4 fm. The scission point of the Th 
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Figure 3: Dissipated energy as function of the internuclear distance R for four 
tunneling velocities. Full line dR/dt = 1 x 104 fm/fs, dashed line dR/dt = 1 x 10° 
fm/fs, dot-dashed line dR/dt = 1x 10° fm/fs, and dotted line dR/dt = 1 x 10" fm/fs. 
Figure adapted from reference [159] 
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Figure 4: Dissipated energy at scission as function of the tunneling velocity v = 
dR/dt. Dot-dashed line neutron dissipated energy, dashed line proton dissipated 
energy and full line total dissipated energy. Figure adapted from reference [159] 
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nucleus is located somewhere at R19 fm. Details about these calculations can be 
found in Ref. [149]. 


The time-dependent pairing equation were solved beginning from the ground 
state up to the scission configuration spanning an interval of tunneling velocities 
v = dR/dt comprised between 1 x 10* and 1 x 10° fm/fs. The initial values of the 
single particle densities and of the pairing moment components are taken as the 
BCS amplitudes for the stationary state. The variation of the dissipated energies 
are plotted in Figs. 3 and 4 for a pairing model with a constant strength G. First 
of all, it can be noticed that the dissipated energy increases especially in the region 
of the outer fission barrier, between R=14 and R=18 fm. This increase of the dis- 
sipated energy in the region of the second barrier is more pronounced in the case of 
state dependent pairing, as it can be seen in Ref. [159], but the behavior is similar 
to the case presented in this work. As expected, in both formalism the dissipated 
energy remains unchanged after the scission point. After scission the rearrangement 
of the single particle levels ceases, and the time dependent pairing equations be- 
come stationary. Therefore, the dissipated energy remains constant. Because, many 
values of the dissipated energy are tested, it is now possible to select the collective 
velocity that pertain to the dissipated energy which infers the experimental data. 
The dissipated energy for the 1x10* fm/fs tunneling velocity is negligible. That is, 
for a very slow deformation of the nuclear system, the nuclear fission process be- 
haves adiabatically. A such value of the dissipated energy corresponds to cold fission 
processes, in which the excitation energy at scission is so small that no neutron are 
emitted. very large. That is, the nuclear systems proceeds through adiabatic states. 
For nuclear velocities that amount to the order of 1 x 10° fm/fs, the dissipated en- 
ergy is very large. That is, the nuclear systems proceeds through diabatic states. 
The dissipated energy amounts to 14 MeV, this values being sufficient to emit at 
least two neutrons from the evaporation from the fission fragments. For tunneling 
velocities larger than 1 x 10" fm/fs, the final dissipated energy decreases and arrives 
on a plateau. This happens because we reached the limit of validity of our model. 
For these large velocities, the passage from the ground state to the scission point is 
too rapid, the occupation probabilities cannot be modified and the system arrives 
at scission with occupation probabilities close to the initial ones. Statistically, the 
Fermi energy of the a nuclear system is roughly about 40 MeV. Accordingly, an 
average nucleon velocity of 8x10’ fm/fs can be associated to the Fermi energy of 
the nucleons. To allow a rearrangement of the nuclear states, the internuclear ve- 
locity should be much smaller than the velocity of the nucleons inside the nucleus. 
Otherwise the nucleons have not the possibility to accommodate with the variations 
of the mean field. 


In Ref. [149] it is shown that for the two pairing mechanisms investigated, the 
dissipated energy have a maximal value around a velocity of 1 x 10° fm/fs. The 
striking difference between the two models is the fact that the dissipated energy 
amounts up to 70 MeV in the case of the state-dependent pairing mechanism. 
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As mentioned above, the cold fission is a process in which the dissipated energy 
is so low that no neutrons can be emitted. So, a dissipated energy of 6 MeV, ap- 
proaching the binding energy of the neutron, can be considered as an upper value 
which characterizes the cold fission processes. {From the Fig. 4 where the depen- 
dence of the dissipated energy at scission is displayed versus the collective velocity, 
the corresponding velocities for an excitation energy of 6 MeV can be inferred as 
1.3x10° fm/fs. It should be noted that is was remarked in Ref. [149] that the 
tunneling velocity is much lower for the state dependent pairing. In fission, the 
intrinsic excitation energies at scission are measured from the number of evaporated 
neutrons from the fission fragments, an it is considered to be around 10 MeV for 
most of the parent nuclei [164]. In this work it is considered that in thermal fission 
the mean dissipation amounts about 12 MeV. For this reasonable value of the dissi- 
pated energy, one deduced that the corresponding velocities range in a small interval 
[4 x 10° - 1.5 x 10° ] fm/fs for a constant pairing strength. A smallest value was 
aso estimated of about 1.2 x 10° fm/fs for a state dependent pairing interactions. 
From the behavior of the deformation energy versus the elongation, a shift of the 
elongation of about 3 fm between the top of the external barrier and the scission 
point was considered. Once this distance available, it is possible to evaluate the time 
for the descent of the barrier. This time has the following values. In order to obtain 
6 MeV dissipated energy at scission, it is T = 2 x 107° s, while in order to have 12 
MeV at scission, it is about 2 x 1077! s. 


We concluded that the tunneling times that can be obtained from calculations 
for the fission processes are strongly model dependent. A realistic pairing interaction 
produce a fission time considerably longer than a constant pairing interaction. The 
fission process could very fast or very slow in superfluid systems to reproduce the 
same value of an observable, depending on the behavior of the pairing interaction. 


Another example concerning the influence of the dynamics of fission is related 
to the pair breaking mechanism. The system of differential equations (1.40) and 
(1.41) can be modified by taking into account a pair breaking mechanism, as done 
in Ref. [157]. In general it is accredited that in cold fission, the distributions of 
the masses of the fission fragments are even-even, due to the the pairing mechanism 
influence [165, 166]. Accordingly, at very low excitation energies, the fragments 
would be fully paired. Very low excitation energies of the fragments are translated 
in high kinetic energies that approach the Q-value of the reaction. A pioneering 
experiment, obtained after improving the experimental procedure, contradicted this 
commonly accepted behavior. In the data of Ref. [167] a first step was achieved: 
no strong even-odd effect was evidenced in the thermal neutron induced fission for 
high values of the total kinetic energy. Surprisingly, even and odd partitions were 
observed experimentally close to their respective Q values for four systems investi- 
gated: 793:235U (nin ,f), 799Pu(nin,f) and 7°2Cf(sf). An even more strange behavior 
was observed in Refs, [168, 169, 170]. The experimental data exhibited higher odd- 
odd fission yields at excitation energies smaller than 4-6 MeV than even-even ones. 
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Figure 5: The dependence of the even-even fission yield as function of the dissipated 
energy is plotted with a full line. The dependence of the odd-odd fission yield as 
function of the dissipated energy is plotted with a broken line. Figure adapted from 
reference [157] 


That lead us to consider the importance of the Landau-Zener effect in the cold- 
fission fragmentation, as mentioned in Ref. [157]. In this case, as a result of the 
interaction in an avoided levels crossing region it is possible to break a pair by 
using dynamical considerations. However, without introducing the dynamics, the 
odd-even structure in fission is explained usually within statistical arguments, as 
for example in Refs. [171, 172, 173]. It is also possible to take into account some 
arguments linked the odd-even structure and to the excitation energy during the 
fission process [175, 176, 177, 178]. 

The fission yields as function of the total excitation energy for the fragmentation 
Kr + '44Ba (even-even), 9?°Rb + !44Cs (odd-odd) is investigated theoretically by 
solving the time dependent pairing equations, generalized by including the Landau- 
Zener effect and the Coriolis coupling. For this reaction, experimental data are 
available for the reaction 73?U(nip,f) [168], allowing us to test the model. 

When the odd-odd partition are treated, the barrier must be larger by adding the 
value of the excitation energy of the unpaired nucleons. That leads to a decrease 
of the fission penetrability. Wheeler introduced this specialization energy in Ref. 
[179]. The fission path was obtained in the framework of the least action principle, 
obtaining the most probable fragmentation. By calculating the single particle levels 
diagrams, the ingredients required to solve the equations of motion are obtained. 
A pairing active level space of 58 levels around the Fermi energy was selected. 
The avoided crossing levels regions where identified. After that, one selected 32 
seniority two configurations for protons and 31 configurations for neutrons that are 
coupled to the seniority zero configuration through avoided levels crossings regions. 
Several values of the internuclear velocity v = R are tested in order to solve the 
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generalized time dependent pairing equations: 5 x 107, 8 x 107, 103, 3 x 10°, 104, 
3 x 104, 10°, 3 x 10°, 10°, and 3 x 10° fm/fs. The ground state configurations of 
the stationary values of the BCS amplitudes were taken as initial conditions to solve 
the equations of motion. The probability to retrieve the system in the ground state 
in a seniority-O configuration was taken as Pop=1. A typical behavior for all the 
seniority-2 configurations was remarked. The general rule exhibits a proportionality 
between the dissipated energy at scission and the collective velocity. 

The final probabilities of realization of all seniority configurations for neutrons 
and protons as function of the internuclear velocities are displayed were obtained 
after solvin the equations of motion. The results evidence a curious behavior: the 
probability of realization of the adiabatic seniority-0 states are close to zero for small 
velocities, exactly in the energy region in which the dissipation is smaller. ; From 
these theoretical results it is possible to understand the main features concerning 
the relation between the final excitation energy and the probability of realization 
of a given channel. If the probability of realization of a seniority-0 configuration, 
that, for the even-even partition is very small at low excitation energies, implicitly 
the probabilities of realization of odd-odd partitions or seniority-2 configurations is 
much larger. The final results are plotted in Fig. 5 by folding on small interval of 
the excitation energy the probabilities of the odd-odd configurations. So, a strange 
experimental result was obtained: at low excitation energy the odd-odd yields sur- 
pass the even-even ones. The even-even yields become larger for excitation energies 
larger that 3-4 MeV. This strange experimental experimental results was explained 
from dynamics. 

Through a dynamical analysis and solving the microscopic equations of motion 
for a fissioning even-even system it is found that the probability to obtain an odd- 
odd partition overcomes the probability of an even-even one at excitation energies 
smaller than 4 MeV, for the same division in mass numbers. The theoretical results 
are in accordance with the experimental behavior of the odd-even distributions at 
high kinetic energies. It is the first time that this behavior was explained within a 
quantum-mechanical approach. 


5 Instead of conclusion 


Today, an increased interest is focused in developing new cycles for nuclear power 
plants which is determined especially by social needs linked to security. It is crucial 
to ensure the safety of nuclear power plants and to handle with responsibility the 
nuclear wastes. As mentioned before, the nuclear fission is still unsatisfactory under- 
stand. A better knowledge of the fission mechanism should be realized in order to 
reduce the risks of the nuclear power plants and to design the new generations. One 
solution is to improve our knowledge on the Th fission. Although thorium itself can- 
not support a nuclear chain reaction, by bombarding thorium to a flux of neutrons 
inside a nuclear reactor, converts this element to uranium 233, which can support 
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fission. For this reasons, the design of nuclear plants have long time considered the 
possibility of a fuel given by a combination of thorium with a fissionable isotope, 
which would initiate the reaction. For this purpose, a better modeling and exper- 
imental determination of cross-section evaluations were needed in the actinide and 
sub-actinide region. As mentioned, actual evaluations are mainly phenomenological, 
the heights of the double fission barrier being determined empirically in accordance 
with a given parameterization of the nuclear level density. This procedure does not 
include the resonance due to the beta vibrational and rotational coupling adequately. 
Also, fission of actinides as well as sub-actinides is not yet understood sufficiently for 
increasing energies above a few MeV. Some efforts must be also involved in a better 
experimental determination of nuclear data by improving the experimental concepts 
and configurations. Fundamental research in the field of nuclear physics is necessary 
to improve our knowledge. The final goal is to obtain values of the neutron-induced 
cross sections within an accuracy of several percents while actually the error bars are 
of almost 30%. So, any theoretical project in fission dynamics is according to such 
efforts. In principle, apart a developments of the theoretical models, the nuclear 
research is also focused towards a better determination of nuclear properties of nu- 
clei by understanding better the nuclear fission mechanism in order to provide more 
reliable data. The economic importance for fission can be evidenced briefly in the 
following. The production of nuclear energy requires a combustible able to furnish 
energy through fission. Between the nuclei constituting the combustible we must 
distinguish among two main types: the fissile nuclei able to fission after a thermal 
neutron capture and the fertile nuclei that lead to a fissile nucleus after the capture 
of a neutron followed by several beta disintegration. The main fissile nuclei are 2°7U, 
235), 239Py and 74!Pu. The main fertile nuclei are 78?Th and 788U. The last ones 
are abundant in the nature and the extraction is not very difficult. Among the fissile 
nuclei, only ??°U can be found in nature. Actually, the U-cycle is imposed as the 
most important way to produce nuclear energy due to historical consideration (the 
first nuclear reactor constructed under the guidance of Enrico Fermi in 1942) and 
due to the development of nuclear weapons. The combustible is obtained by mixing 
two isotopes of uranium (235 and 238). Some drawbacks are raised by the U-cycle. 
The first one is linked to the criticality of the actual reactors. The consequences 
of the Tchernobyl accident were proved disastrous. Secondly, the administration of 
the nuclear wastes is another important problem. Thirdly, the resources are limited, 
the abundance of ??°U is only 0.72% of the natural uranium Improved reactor types 
are needed because uranium resources are available only for 30 y with the actual 
utilization [180]. In such conditions, another cycle based on a mixing of 7°7U and 
232Th becomes recently very promising for the nuclear energy production, namely 
the Th-cycle. Unfortunately, this cycle has also some drawbacks. A lot of invest- 
ment and fundamental researches must be realized in the field. The project of such 
hybrid reactor is launched at European and International levels. As mentioned in 
the Ref. [181] a complex European project started to propose a technological route 
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to reduce the risks associated to nuclear waste, route based on transmutation in 
ADS. it is necessary to determine accurately the fission and capture cross sections 
for 793Pa, 78°Th, ?32Th and the trans-uranium elements. These cross sections must 
be known at least within an accuracy of 15%. Actually, the evaluations obtained 
for these nuclei give a precision smaller than 30%. In the case of ?°?Th the cross 
section accuracy is required to be about of 2 % and in the case of 7°°U a value 
of 1% is necessary. This last isotope supplies the main exothermic reaction of the 
cycle. More over, the actinide data are scarce and often discrepant, especially for 
232Th and 733U and must be supplemented by further experimental studies. Due 
to the lack of knowledge in the field, most countries abandoned this Th scheme 
decades ago, but increasing concerns about the diversion of plutonium from spent 
nuclear fuel to the construction of nuclear weapons has prompted reexamination. 
IAEA has maintained an interest in the thorium fuel cycle that offers the previ- 
ous advantages. Momentary, the following fission cross sections must be improved: 
232— 23477 230—234Dq and 230-233Th. The most important requirements can be clas- 
sified as follows: determination of the basic cross-sections, evaluation of available 
experimental data to compute cross-sections and improvements of actual models. 
The largest activities in Europe concerning these goals are coordinated within the 
n_TOF project and cover experimental measurements and cross-sections and model 
evaluations. The activities plan to improve the understanding of the mechanism of 
fission in order to improve the cross section evaluations. 
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